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Abstract 

We have observed the shadowing of galactic cosmic ray flux in the direction of the 
moon, the so-called moon shadow, using the Tibet-Ill air shower array operating at 
Yangbajing (4300 m a.s.l.) in Tibet since 1999. Almost all cosmic rays are positively 
charged; for that reason, they are bent by the geomagnetic field, thereby shifting the 
moon shadow westward. The cosmic rays will also produce an additional shadow 
in the eastward direction of the moon if cosmic rays contain negatively charged 
particles, such as antiprotons, with some fraction. We selected 1.5 x 10^° air shower 
events with energy beyond about 3 TeV from the dataset observed by the Tibet-Ill 
air shower array and detected the moon shadow at ~ 40(7 level. The center of the 
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moon was detected in the direction away from the apparent center of the moon 
by 0.23° to the west. Based on these data and a full Monte Carlo simulation, we 
searched for the existence of the shadow produced by antiprotons at the multi- 
TeV energy region. No evidence of the existence of antiprotons was found in this 
energy region. We obtained the 90% confidence level upper limit of the flux ratio of 
antiprotons to protons as 7% at multi-TeV energies. 

Key words: Cosmic rays, Air shower, Anti matter. Observation and data analysis, 
Simulations 

PACS: 96.50.sd, 94.20.wq, 25.43.+t, 95.30.-k, 95.75.-z, 98.80.-k 



1 Introduction 

Observation of antiproton abundance in cosmic ray flux raises the possibility of 
a baryon-symmetric universe and the propagation of cosmic rays in interstellar 
space. Recent measurements of antiproton flux, however, appear to be almost 
within the conventional cosmic ray physics in which antiprotons are produced 
as secondary particles of cosmic ray interactions with interstellar gas. For 
example, antiprotons are produced mainly by collisions of cosmic-ray protons 
with interstellar hydrogen gas as p+p — > p+p+p+p. Accelerator experiments 
measured the production rate (p/p) of antiprotons to protons to be on the 
order of 10~^ at energies greater than 10 GeV. 

Recently, measurements of absolute flux of antiprotons below a few GeV were 
carried out using the magnet spectrometer equipped with a track detector to 
identify the charge and momentum of each incident particle [1,2,3,4]. Among 
these, the CAPR1CE2 experiment extended the energy range of the spectrum 
up to about 50 GeV. It is obvious that accurate measurements of p flux are 
key to testing current propagation models of cosmic rays in the Galaxy. Of 
course, a p "excess" from the reliable propagation model may lead to discover 
possible sources of primary antiprotons such as dark matter annihilation and 
evaporation of primordial black holes. Strong ct al. [5] made a detailed calcu- 
lation of the antiproton flux, diffuse gamma rays and other cosmic ray fluxes 
and compared it with the recent BESS results available in the energy region 
below a few GeV [4] . They found that the conventional local cosmic ray mea- 
surements, simple energy dependence of the diffusion coefficient, and uniform 
cosmic ray source spectra through the Galaxy fail to reproduce simultaneously 
both the secondary to primary nuclei ratio and antiproton ffux in this energy 
region. 

* Corresponding author. 

Email address: hibino@icrr.u-tokyo.ac.jp (K. Hibino ). 
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The experiments mentioned above also present the flux ratio, p/p, of antipro- 
tons to protons in the cosmic rays. It seems that the observed results are within 
the scope of the prediction of the standard leaky box model [6]. Some results 
[1,3], however, show a tendency that the p/p ratio increases with increasing 
primary energy, although the amount of data is insufficient. 

In the simple leaky box model [7], the ratio of antiproton flux fp{Ep) to proton 
flux fp{Ep) is calculated as [8] 



where Agsc = pl^CTesc (c is the light velocity) is the mean amount of matter 
(density p) traversed by a nucleon of velocity /3c, Tesc is the mean time spent 
by the cosmic rays in the confinement space, and A at is the mean free path of 
the nucleon in the interstellar matter. In addition, Zj^^ is the production rate 
of antinucleons by nucleon - nucleon (hydrogen atom) interactions, which is 
almost constant when Feynman scaling holds, while weakly depending on the 
spectral index of cosmic rays. According to the HEAO-3 data on the B/C ratio 
in the primary cosmic rays up to about 50 GeV/n, X^sc oc and 5 ~ 0.6 [9]. 
Because the mean free path of the nucleon, Aat, in this energy region is almost 
constant, the flux ratio oip/p should decrease with increasing primary energy. 
At high energies greater than 100 GeV, however, a value oi S ^ 0.3 is expected 
because the cosmic ray anisotropy is less than 10^^ in the energy region below 
100 TeV. Actually, 6 should take a value of 1/3 when the magnetic held in 
interstellar space has a Kolmogorov-type spectrum. 

A closed galaxy model [10] predicts a considerable enhancement of secondary 
antiprotons. In this model, the cosmic ray flux observed locally in the solar 
system comprises young and old components. Cosmic rays in the young com- 
ponent show almost identical behavior to that discussed above, whereas the 
old component is contained in the galaxy halo. The old component of cos- 
mic rays, however, consists only of protons because heavier nuclei in the old 
component are completely broken up during their long conflnement time in 
the halo. Consequently, the value of Agsc for protons in the old component 
becomes very large. This model can engender a considerable enhancement of 
secondary antiprotons because of the existence of such old protons [11]. 

Until now, almost all measurements of antiprotons are limited in the energy 
region below a few GeV because of its extremely low flux and weight limit 
of the instrument on board the balloon. In this energy region, the interpre- 
tation of observed spectra is not free from uncertainties about the secondary 
antiproton production because of the nuclear threshold effect and because of 
uncertainties in the solar modulation effect. For those reasons, it is still dif- 
flcult to distinguish which model is the best for explaining the antiproton 
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spcctrTim. Considering the latest research, the possibihty of the existence of 
antiprotons in the TeV energy region is shght, but we cannot exchide an exotic 
model such as an anti-galaxy production model yet (see the review paper by 
Stephen & Golden [12]). Therefore, it will be important to measure p/p ratio 
in the TeV energy region. The results of such measurements might enable us 
to investigate the possible production of antiprotons from dark matter at high 
energies [13]. 

Urban et al. [14] first discussed the possibility of observing antiprotons using 
the moon shadow; the first observation was done by the Tibet-I experiment 
[15] by observing the sun shadow, giving an upper limit of the p/p ratio at 10 
TeV [16]. The Tibet air shower experiment has been observing the shadowing 
of galactic cosmic rays continuously by the moon and sun since 1990 [17]. We 
began construction of the Tibet-Ill air shower array in 1999 by increasing the 
number of detectors; we completed it in Fall of 2003. The present Tibet-Ill 
air shower array consists of 789 scintillation counters of 0.5 m^, each of which 
is placed in a 7.5 m square grid, covering a total area of about 37000 m^ [18]. 

In this paper, we present the result on the search for antiprotons based on the 
moon shadow observed during the period from 1999 through 2004 using the 
Tibet-Ill air shower array. 



2 Experiment 

The Tibet-Ill air shower array used in this experiment was constructed in 1999 
at Yangbajing (4300 m a.s.l.) in Tibet and operated until 2004. The array, 
corresponding to the inner part of the full-scale Tibet-Ill air shower array, 
consists of 533 scintillation counters covering 22050 m^ [19,20]. The mode 
energy of detected events is about 3 TeV for proton-induced showers and the 
angular resolution is 0.9°. The systematic error of the energy determination 
of primary particles and systematic pointing error of the array have been 
well calibrated by comparing the observed displacement of the moon shadow 
because of the geomagnetic field with the Monte Carlo simulation, as discussed 
in the paper [21]. 

With the advent of this array, we observed 6.1 x 10^° events during the period 
from November 1999 through December 2004 (1041 live days in total). These 
events were selected by imposing the following conditions: 1) Each shower 
must fire four or more counters recording 1.25 or more particles; 2) all of fired 
counters or eight of nine fired counter which recorded the highest particle 
density must be inside the fiducial area; and 3) the zenith angle of the arrival 
direction must be less than 40°. After the these selections, 1.5 x 10^° events 
remained for further analyses. 
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3 Simulation 



We have done a detailed Monte Carlo simulation to estimate the flux of an- 
tiprotons from the shape of the moon shadow by cosmic rays under the in- 
fluence of geomagnetic field. A deviation of the observed shadow from the 
simulated one based on the normal primary cosmic rays (without antipar- 
ticles) would give the proof of the existence of antiprotons pouring on the 
Earth. For the normal primary cosmic rays, we used a best-flt curve for each 
component based on the experimental data measured by direct observations 
on board the balloon and satellite (for experimental data, see the paper by 
Gaisser et al. [22] and references therein). 

We used the CORSIKA Ver. 6.200 code [23] with the QGSJET interaction 
model for generation of air showers in the atmosphere. All secondary particles 
produced by a primary particle were traced until their energies becomes 1 
MeV in the atmosphere. Simulated air-shower events were then input to the 
detector with the same detector configuration as the Tibet-Ill air shower array 
using the Epics uvS.OO code [24] to calculate the energy deposit of these shower 
particles and were stored in the disk with the same format as the experimental 
data. 

For the geomagnetic field, we used the Virtual Dipole Moment (VDM) model 
[25] in which the dipole moment at altitude > 600 km is taken to be 8.07 x 
10^^ G • cm^. We also made the simulation based on the International Geo- 
magnetic Reference Field (IGRF) model [26] to compare with that based on 
the VDM model. The difference of the defiection of particle trajectries between 
the two models is then estimated to be within 5 ~ 7% in the part of the sky 
where the moon is visible by the Tibet-Ill air shower array. We found, how- 
ever, that the IGRF model gives almost the same result as the VDM model 
for the actual moon shadow simulation in which the angular resolution of the 
array is taken into account. Hence, we adopted the results obtained based on 
the VDM model in this paper. 

The simulation for the moon shadow was performed through the following 
steps: 

(1) Primary particles in the energy range from 0.3 TeV to 1000 TeV are 
thrown isotropically toward the observation site on the top of the atmo- 
sphere along the revolution orbit of the moon. Air shower events gener- 
ated by these particles in the atmosphere are then collected within the 
circle of 300 m radius from the center of an air shower array. This array 
has the same structure and detector response as the Tibet-Ill air shower 
array. 

(2) All air shower events triggered by the air shower array are analyzed in the 
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same way as the experimental data to estimate the size and the arrival 
direction of each event. 

(3) An opposite charge is then assigned to the primary particle of each ana- 
lyzed event; and these particles (anti-particles) are again shot back uni- 
formly in the ±10° x ±10° direction centering on the moon: they origi- 
nated from the air shower array at Yangbajing and tracked through the 
geomagnetic field up to the moon. 

(4) We select only the events hitting the moon; these selected events are again 
traced back from the moon to the Earth after reversing the charge of each 
primary particle. The final direction of each cosmic ray is determined by 
smearing the difference between the moon direction and the estimated 
air shower direction where the angular resolution of the array is taken 
into account. 

Using this method, we can simulate the moon shadow that is capable of direct 
comparison with the observed data. 



4 Data Analysis 



For analysis of the moon shadow, an on-source event is defined as an event 
which falls on the ±5° x ±5° square window centered at the moon direction. 
The background is estimated using the events coming from the following eight 
off-source windows of the same size as on-source. These positions are located 
on the same zenith angle as the moon, but separated by ±5°, ±10°, ±15°, and 
±20° in the azimuth angle. The number of events averaged over these eight 
off-source windows was taken to be the background [17]. 

Shown in Fig. 1 is the moon shadow observed using the events obtained for 
1041 five days with the Tibet-Ill air shower array. The significance at the 
peak of the shadow is estimated as about iOa level and the shadow is shifted 
westward by about 0.23°. The mode energy of observed events is estimated as 3 
TeV for proton-induced showers. Note that a 1 TeV proton is shifted westward 
by 1.6° at Yangbajing. This means that a major effect of antiprotons should 
appear in the eastward angular distance smaller than about 0.5° on the moon 
shadow observed with the Tibet-Ill air shower array. The systematic error of 
the energy estimation of primary particle is estimated to be ± 8% level and 
the systematic pointing error of the array is smaller than 0.011° [21]. Using 
these basic data and simulation results, we discuss the flux of antiprotons in 
the next section. 



7 




Angular Distance (degree) 



Fig. 1. Significance map of the deficit event densities observed with the Tibet-Ill air 
shower array for 1041 live days. The origin of the coordinate is taken in the direction 
of the moon. The white dashed circle denotes the apparent size and position of the 
moon. A right-hand side scale expresses the level of significance of the deficit event 
density in terms of the standard deviation a. The significance at the central peak 
of the shadow is about 40a on the two-dimensional analysis [17]. 

5 Results and Discussion 

Figure 2 shows the one-dimensional distribution of the observed deficit events 
around the moon in the north-south direction, where the vertical axis rep- 
resents the number of deficit events inside the range of ±1° in the east-west 
direction. The origin of the coordinate is taken to indicate the apparent direc- 
tion of the moon. The position of the shadow center is almost consistent with 
that of the moon because the east-west component of the geomagnetic field 
is negligibly small at Yangbajing latitude, so that the cosmic rays are not de- 
flected in the north-south direction. As this figure shows, our data show good 
agreement with the MC simulation {x^/d.o.f. ~ 1.15). The peak position of 
the distribution in the north-south direction is estimated as 0.0048° ± 0.011°. 

Figure 3 shows the one-dimensional distribution of the observed deficit events 
around the moon in the east-west direction, where the vertical axis represents 
the number of deficit events inside the range of ±1° in the north-south di- 
rection. The observed moon shadow is shifted westward by 0.23°, which is 
consistent with the MC simulation. There exists, however, no evidence indi- 
cating deficits of cosmic rays at the opposite position around 9 — 0.23° in 
the eastward direction, corresponding to the particles with negative charge 
(anti-matter) such as p. He, C,.., and Fe, if any. 

Next, to evaluate the contribution of antiprotons to the moon shadow, we 
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Fig. 2. One-dimensional distribution of the observed deficit events around the center 
of the moon in the north-south direction. The vertical axis represents the number 
of deficit events inside the range of ±1° in the east-west direction. The open circles 
denote the experimental data and the filled triangles denote the MC data. The 
solid line represents a best fit curve according to the function fi{9). Each error 
bar denotes ±.la. The histogram of the following figure shows the residual of the 
observed deficit event distribution to the function fi{0). 

obtained two kinds of moon shadow that were cast by all cosmic rays and 

protons, respectively. We assume that the deficit event distribution in each 
shadow can be expressed as a superposition of Gaussian functions. It is found 
that four Gaussian functions are adequate for fitting to both distributions in 
the angular distance smaller than about 5° from the center of the moon. 

The moon shadow by all cosmic rays is then expressed as 



f^n, 4-, I 41n2x(g-M„,,,)^ >| 

/l(^) = 2^^aH,iexp 2 . 

1=1 \ '-^all,i J 



The moon shadow by protons is also expressed as 



hie) = Y: a,, exp '-^ ^ , (3) 

1=1 \ ^p,i ) 

where Q is the angular distance from the moon center in the west-east direc- 
tion and M and a are the fitting parameters for the distribution function, 
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Fig. 3. One-dimensional distribution of the observed deficit events around the center 
of the moon in the east-west direction. The vertical axis represents the number of 
deficit events inside the range of ±1° in the north-south direction. The open circles 
denote the experimental data and the filled triangles denote the MC data without 
antiprotons. The solid line is a best-fit curve according to the function f3{9). Each 
error bar denotes ztla. The histogram of the following figure shows the residual of 
the observed deficit event distribution to the function f3{0). 

respectively, for all cosmic rays and protons. 

On the other hand, the observed moon shadow should be expressed by the 
following function /3(^) 

f3{e) = am + bM-e), (4) 



where the first term represents the deficit in cosmic rays and the second term 
represents the deficit in antiprotons. From our simulation, the fraction of pro- 
tons to all cosmic rays is estimated as 62 ± 1% of all cosmic rays, so the ratio 
6/0. 62a corresponds to the p/p ratio. As shown in Fig. 3, the observed deficit 
events can be fitted by the function fsiO) with the parameters of a = 1.64±0.10 
and b = —0.10 ± 0.08 {x^/d.o.f. ~ 1.50), where the fitting was made under 
the boundary condition that the integral of afi (9) should be smaller than the 
total number of observed events. 

To calculate an upper limit of the fraction of antiprotons, the confidence in- 
terval of the parameter b is estimated using Feldman and Cousins Statistics 
[27] because parameter b is negative. The upper limit of the p/p ratio at the 
90% confidence level was then obtained as 0.07. The present result is plotted 



10 




I • A 



V...-- 



Present work 



H Buffington (1981) 

O Golden (1984) 

O Bogomov (1987-1990) 

* Streitmatter (1990) 
^ Salmon (1990) 

• IMAX(1996) 

□ CAPRICE (1997) 

A Basin! (1999) 

■ BESS (1999) 

O BESS (2000) 

H CAPRICE2 (2001) 

* HEAT (2001) 

V MACRO (2003) 

T L3+C(2005) 

O TIBET-I (1991-1993) 

• TIBET-III (1999-2003) 



102 



103 



lO"* 



Primary Energy (GeV) 



Fig. 4. The antiproton/proton ratio at the top of the atmosphere. The experimen- 
tal results are as follows: Buffington (1981) [28], Golden (1984) [29], Bogomolov 
(1987-1990) [30], Salomon (1990) [31], IMAX [32], CAPRICE (1997) [33], Basini 
(1999) [1], BESS (1999) and BESS (2000) [4], CAPRICE2 (2001) [3], HEAT (2001) 
[2], L3+C (2005) [34], Streitmatter (1990) [35], MACRO (2003) [36], L3+C (2005) 
[34], and opened circles are those by TIBET-I (1990-1993) [16]. The present result 
by TIBET-III (1999-2004) is denoted by the filled circle. The dashed lines, upper 
and lower limit, are the calculation based on the leaky box model by Simon, Molnar 
&: Roesler [6]. The dotted lines are for a model that includes both leaky box and 
the presence of antistars. The rigidity dependent confinement of cosmic rays in the 
Galaxy is assumed to be oc R^^, and two curves A and B are the cases of 5 = 0.7 
and 0.6, respectively [12]. 



in Fig. 4, together with other results. 



In this figure, four experimental data are plotted at multi-TeV energies. Among 
them, two are from the Tibet air-shower experiment using the moon shadow 
(present work) and the sun shadow (Tibet-I [16]). The present result gives the 
most stringent constraint on the p/p ratio compared to other data. 

In conclusion, we showed that the moon shadow observation can provide 
unique data on the p/p ratio at the multi-TeV region. A further observa- 
tion of the moon shadow and a fine tuning of the MC simulation will give a 
more stringent constraint to the flux ratio of antiprotons to protons in the very 
near future. If these data are combined with the new data that are obtain- 
able in the energy region around 100 GeV by the PAMELA satellite mission 
(launched on board the Russian satellite in 2006) [37] and by the forthcoming 
space experiment AMS [38], we might be able to get new information about 
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the production of antiprotons in the TeV energy region. 
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